
function [tt , xx] = Cr3bp_Orbit2D(tVec, x0, aux)
%
% 画CR3BP轨道，2D
%
% 时间：2024年06月03日
% 作者：张晨
% 单位：CSU
% 邮箱：chenzhang@csu.ac.cn
%%%%%%%%%%%%%%%%%%%%%%%

mu = aux.EMRot.mu;

options = odeset('Reltol', 1e-8, 'AbsTol', 1e-8) ;
[tt , xx] = ode113(@Cr3bp_Eqm2D , [tVec(1) , tVec(2)] , x0 , options , mu);

end
